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一 类 三 次 有 理 插值 样 条 的 逼近 性 质 * 


刘 琳 ， 唐 月 红 
(南京 航空 航天 大 学 理学 院 ， 南 京 210016) 


摘 要 : 曲线 的 保 形 播 值 是 几何 外 形 设 计 的 重要 课题 。 本 文 构造 了 一 类 带 控制 参数 且 包 含 极点 
的 (3,2)* (k = 1,2) 阶 有 理 插值 样 条 。 对 于 给 定 的 单调 和 保 凸 数组 ， 通 过 对 样 条 中 参数 的 适 
当选 取 达 到 保 形 的 旭 的 。 对 于 (3,2)* (k = 1,2) 阶 插值 曲线 的 形状 控制 问题 进行 了 研究 ， 推 
导出 了 将 此 插值 曲线 约束 在 给 定 的 折线 和 二 次 曲线 之 上 、 之 下 或 之 间 的 充分 条 件 。 最 后 本 文 
以 Peano-Kernel 定理 为 工具 ， 讨 论 了 该 插值 的 逼近 性 质 。 给 出 的 数值 例子 说 明 这 些 方法 的 有 效 
性 。 
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对 于 [a, 四 上 的 一 个 分 划 A :a = xo € o < zn = bH (ru D WEG BURE A ERU PROS 
A Rio MXR. EK, SKAER (3, 1)*, (3,2)*(k = 1,2) 的 研究 在 理论 和 应 用 上 已 
取得 了 一 些 结果 [23， 但 包含 极点 的 有 理 方法 很 少见 。 文 献 四 给 出 带 极点 的 有 理 插值 ， 但 不 是 
样 条 。 本 文 研究 一 类 三 次 有 理 样 条 (3,2)* (k = 1,2)， 即 满足 下 列 条 件 的 有 理 样 条 函数 R(z) 


e? (1-8)? 


'lrg-09) 1440 S z € [wi-1, zi], 


R(x) = span{1, 0 


R(z) € RÜNA), R(z)e RÓ(A), ze [a,b], 
R(zi) = fi 1—0,1,::-,mn, 


其 中 f; 为 型 值 ，g; > 0 为 控制 参数 ，0 = (r— zi; 1)/hi, hi = zi — zii 

用 多 项 式 研 究 保 形 插值 是 比较 困难 的 ， 如 用 分 段 三 次 多 项 式 构造 C1 保 凸 插值 条 件 就 十 分 其 
刻 ， 即 使 对 一 般 的 凸 数组 也 需要 插入 较 多 的 内 节点 才能 满足 条 件 。 样 条 插值 是 曲线 和 曲面 设计 
中 的 有 力 工具 。 一 些 作者 已 经 研究 了 不 少 类 型 的 样 条 函数 用 于 曲线 和 曲面 的 形状 控制 。 近 些 年 
来 ， 有 理 样 条 ， 特 别 是 有 理 三 次 样 条 以 及 它们 在 形状 控制 中 的 应 用 已 有 较 多 研究 。 但 在 播 值 条 
件 固 定 的 情况 下 ， 普 通 的 插值 样 条 曲线 如 三 次 插值 样 条 ，B- 样 条 等 ， 因 为 插值 格式 固定 而 无 法 
控制 曲线 的 形状 。 本 文 基于 构造 的 有 理 插值 函数 ， 研 究 保单 调和 保 凸 有 理 样 条 播 值 方法 ， 并 讨 
论 了 这 种 插值 曲线 控制 于 两 给 定 曲线 的 问题 以 及 它 的 逼近 性 质 。 
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一 类 三 次 有 理 插值 样 条 的 逼近 性 质 


定理 2.1 (3,1)! Wr R(z) 可 以 表示 为 


R(z) 


其 中 R'(;) = di, i-1lee 


= fiii (0) + fip2(0) + hidi_193(0) + hidia(0), 
sno di 为 B 由 参数 。 


(eLO) oO) (0 aO) = (1.8.5 7) 
Teo r dH wm 
-E Oa cB. Onus 
事实 上 ， 定 理 2.1 的 证 明 是 明显 的 。 我 们 可 直接 验 明 pi(6) 满足 
pili — 1) = pi 1) = fij PG 1) = pil 1) =0, 453-12, 
且 wi(b) € Rs,2。 此 外 ， 容 易 验 证 
pil) = 92(1 — 0), 3(0) = —pa(1 - 6). 
对 于 (3,2)? 阶 有 理 插值 样 条 函数 仍 可 表示 为 (1)， 但 由 于 在 内 节点 上 有 
RE) (zi +0) = R® (z; -0), k=0,1,2, 
这 时 d; 需要 满足 方程 组 
Hiyidi—1 + (o; Hi + o441Hi)di + Hidi41 = ci, 
其 中 
Hi = (o]- h/f Bi= 3 4qi  2q7, 
G= |H: (1+ ais pŠ fg +o A es a 


类 似 地 ，(3,2)2 阶 R(z) 也 可 表示 为 
fiz 193 (0) +t Fio, (0) T h2Di ig3(0) t k? Diga(0), 


(pı (8), v3 (0), P3(0), 7a(9)) 


TE [zi-1, zi], 


R(z) = 


10-5! 0 


2 - y pg 
三 (1, 6, pue s. y 9-9) x 
工 十 qi(1 zd 8) 1+ qi 0 0 0 Bm? 


0 0 
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(2) 


(3) 


(4) 


(6) 
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而 R"(z,) = D, 满足 下 列 方程 组 


HiDi i+ (oH; + o 41H 3)Di + Hii Dii = Ci, (8) 
这 里 
T. S fin - fi fi | hi-fi fii 
H; = hb; , Ci= EUER IER 
由 以 上 我 们 可 以 把 三 次 有 理 插值 样 条 R(z) 表示 为 
= 20), 
其 中 


pi(z) = (1—0)3fi1+0(1— 0) (rifi-1 + hidi-1) + (1 — 08? (rif; — hidi) + 6? fi, 
qí(z) = (1-0? --(1—0)0(r; 1) -?, ri 22 gi (qi 9-1), 
并 有 Ti > 3。 
3 ” 保 形 插值 的 有 理 样 条 方法 
对 于 具有 单调 递增 性 质 的 插值 数据 点 ， 即 有 < fo <… < fus 或 
A= Gi - fi) ue 20, i-12,-.,n- 1, 


有 以 下 定理 。 
定理 3.1 对 于 给 定 的 {(zi, fi di), i=0,1, n} {fi} 为 单调 递增 数组 ，(3,2)? Br R(z) 
在 [zi-1, zi] 上 为 单调 递增 的 充分 条 件 为 


a da, g} < v3, 3 Ai; #0 时 ， 


(10) 
di-1 = di = 0, 当 A; = 0 时 ， 
必要 条 件 为 di > 0, i= 0,1,.… ,n。 
证 明 必要 条 件 的 证 明 是 显然 的 。 只 须 证 (10) 为 充分 条 件 。 由 (9) 求 一 阶 导数 并 化 简 整理 得 
Mi(z) 
q(x)? 


R'(z) = 
其 中 Mi(x) = a0* + b6? + c0? 4- d0 -- e. M 
a = A;(3 — 4r; + r2) — (ri — 3)(di-1 + di), b = 2(r; — 3) (dii + di) - 2Ai(3— 4r; + 72), 
c= A;(3 — Ari +r?) — (ri — 6)(di-1 rj) - 2r;A;, d= 2r;A; — 4di-ı — 2di;, e= di1- 
40-—.,520, I 


As! +B Ds+E 
Mic) sux) ZE DIM But E ENO E ud. 


B =b+2c+3d+ 4e = 2(riA; 一 di-1), C = c + 3d + 6e = (3 + r2)A; — ri(di-1 十 di), 


D =d + 4e = 2(ri A; — di), E =e = dii. 
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因为 R'(z) 的 分 母 大 于 零 ， 则 尽 (z) > 0 等 价 于 Mi(z) > 0， 而 Mi(z) > 0 成 立 只 须 
A20, B20, C20, D20, E20, 
故 三 次 有 理 样 条 插值 函数 R(x) E x € [ei Xi] 上 为 单调 递增 的 充分 条 件 为 
ri; — dii 2 0, 


(3 - 02)A; — ri(di-1 + di) 2 0, 
Tii m di 2 0. 


当 A; 关 0 时 ， 这 组 不 等 式 等 价 于 


因为 7; > 3， 要 使 上 式 成 立 ， 只 须 


当 A; = 0 时 ， 结 论 是 显然 的 。 
EEG 对 于 单调 递减 的 {fi}?-o， 令 == RB. 


考虑 严格 西 的 数据 点 {xi fii = 0,1,- ,n}, BIA, <1 A5 «e An 1 < Ane 为 避免 插 
值 曲线 中 出 现 直 线段 ， 则 要 求 数据 点 满足 


do <Ai<di<:..<di1<Ai;i<di<:…<An < ds. 


定理 3.2 ”对 于 给 定 的 {(zi, fi, di),i = 0,1,… np Uo APERIRE, (3,2)! Br R) E 
[zi-1, 2i] 上 为 凸 的 充分 条 件 为 


rmn E icta) 


i i—1 


(11) 
证 明 将 (9) 求 二 阶 导数 化 简 整理 得 


n,a AS +B? c Cs - D 
POS hatt” 
类 似 于 定理 3.1 的 证 明 ， 当 满足 下 列 不 等 式 组 


(di-1 — di) + ri(di ~ A;) È 0, 
di E Ai 之 0， 
(12) 
—(di-1 — Ai) > 0, 
(di-1-— di) — ri(di-1 — Ai) 2 0, 


A As Bs 4CstD20, AMA R'(z)20.fmm^5xX(2) TERURITA (11). 
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4 ”约束 插值 问题 


设 f(x) 是 被 插 函 数 ， 令 Ra) 是 在 [zo, zw] 上 由 (1) 式 所 定义 的 三 次 有 理 插值 样 条 函数 ，g(zx) 
EE [zo tn] 上 定义 的 线性 函数 ， 或 以 zo < zl < … < En- < zn 为 分 点 的 分 段 线 性 函 
数 ， 数 据 {(zifi, di), = 0,1, n} WÆ fi > (或 <)g(zi), i = 0,1,… ,n。 如 果 对 所 用 E 
[zo, zn] A R(z) > (IÈ <) g(x), i = 0,1,… ,n, WU R(z) 称 作 被 约束 于 g(x) 之 上 (下 ) 的 三 次 有 
理 插值 样 条 函数 ，g(z) 被 称 作 下 (上 ) 约束 曲线 。 

定理 4.1 给 定数 据 {(zi, fi di, gi) i 20,1--,nb Hf; € gi» XXRgi = g(z;), W 
由 (1) 式 定义 的 三 次 有 理 插值 样 条 函数 R(z) 在 [xz;_1,zi] 上 位 于 直线 段 g(z) 之 下 的 充分 条 件 
Æ: 参数 ri (r > 3) 满足 如 下 不 等 式 组 


| rifiit+hidii— gi — gi(ri—1)<o, as 
rif; — hidi — gií(ri - 1) - gi-1z 0. 
证 明 由 (9) 式 容易 看 到 ， 在 [zi ci oi] 上 
a(z)- 0-9? 0-091) 20, R()- PD < gla), 
等 价 于 pi(z) 一 aí(z)g(z) < 0s $ Zi(z) = mi(z) 一 qi(z)g(z)， 于 是 
Zil) = (1-0) fi-1 4 6(1— 0 (rif 1 + hidi-i) + (1 — 0)0? (r: fi — hidi) 
H8 f, — (a — 6)? + (1— 6)6(r; — 1) - €] [1 — 0)gi-1 + 6gi]. (14) 


因为 
[a — 0)? (1 — 6)6(r; — 1) - €?] [(1 — 0)g;-1 + 0g;] 
= (1- 0i + (1 — 9*6[gi + gi-i(ri — 1)) + (1 — 0)0? [giri — 1) + gi-i] + 079; 
故 (14) 式 变 为 
Zi(z) = (1 — 0 (fix ~ gi-1) + (1 - 0?6E, + (1— 0) F; + (fi — gi), 
这 里 
E; = rifi-1 + hidi-1 — gi — gii(ri 1), F; = rifi — hidi — gi(ri — 1) — gi-i- 


wR E; «0, F; € 0, WA Zi(£) < 0 对 所 有 z € [ri 2,] IRA, FERE. 
类 似 于 定理 4.1 的 证 明 ， 对 于 “下 约束 ”的 情况 ， 有 下 述 充分 条 件 的 定理 。 
定理 4.2 给 定数 据 {(zi, fandigi) i = 0,1, ,np Hf; » gt. XH gr = 入 (ci)， 则 
由 (1) 式 定义 的 三 次 有 理 插值 样 条 函数 RR(z) 在 [zi;_1, i] 上 位 于 直线 段 9%*(z) 之 上 的 充分 条 件 
是 : 参数 ri (ri > 3) 满足 如 下 不 等 式 组 
| rifi-1  hidi-i — gt — gii(ri—1)>0, (18) 
rif; — hidi — gi (ri - 1) — gi. 2 0. 
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由 定理 4.1 和 定理 4.2 可 得 下 面 “ 上 下 同时 约束 ”的 定理 。 

定理 4.3 ”给 定数 据 {(zi, fi,di, gogh i=0,1, ,nb» Bgt < fi < 9;。 则 由 (1) 式 定义 的 
三 次 有 理 插 值 样 条 函数 RR(z) 在 [mia zi] 上 位 于 直线 段 g(z) 之 下 且 位 于 直线 段 g*(z) 之 上 的 充 
分 条 件 是 ， 参数 mi(r; > 3) 满足 不 等 式 组 (13), (15). 

设 f(z) 是 被 插 函 数 ， 令 R(z) 是 在 z € [zo,zn] 上 由 (1) 式 定义 的 三 次 有 理 插值 样 条 函 
3, ke) Ær € [zo,zn] 上 定义 的 二 次 函数 ， 或 以 Zo < zl < … < xa < zn 为 分 点 的 分 
段 二 次 函数 ，k;, k DIER k(x) XE z; 处 的 函数 值 和 导数 值 。 仿 上 面 的 证 明 ， 可 得 到 将 插值 曲 
线 约 束 于 给 定 的 二 次 曲线 之 下 、 之 上 或 之 间 的 定理 。 

定理 4.4 给 出 数据 {(zi, fi di kiki), 4 0,1, ,nn}， 且 到 <， 则 由 (1) 式 定义 的 三 次 
有 理 插 值 样 条 函数 R) E [eii zi] 上 位 于 二 次 曲线 段 k(z) 之 下 的 充分 条 件 是 ， 参 数 7; (ri > 
3) 满足 如 下 不 等 式 组 


(1+ ri)(fi-1 — ki-1) + hi(di-1 — k;—1) < 0, 
ri(fi-1 + fi) + hi|dii — di — (ri — 1)ki_1] ^ ki — (2r; — 1)ki-ı < 0, (16) 
(ri + 1)fi — hildi + kj 4) 一 (ri 一 1)ki <0. 
证 明 由 (9) 式 容易 看 到 ， 在 [zi r] E 
qi(z) = (1 —0)? + (a — 6)6(r; — 1) - 8? 7 0, 
k;(z) = (1 — )3k; + (1 — 0)0(2ks_1 + KL Lh) + 85k, 
于 是 l 
nr) 
R(z) e qi(x) < k(z), 
等 价 于 pi(z) 一 qi(z)k(z) < 0, 4 Ui(z) = Di(z) — aí(z)k(z). TÆ 
Ui(r) = (1—0) fi.1 - 6(1 — 0) (rifi-1 + hidi) 


(1 — 6)8*(rif; — hidi) — [1 — 6? + (1 — 6)6(r; — 1) + €?] 


[(1 — 0? k;.1 + (1 — 0)0(2k; 1 + kj hi) + 9*k;], (17) 
而 
[(1 ~ 0)?+(1—0)0(r;—1)+0°][(1—0)2k1+(1— 0)0(2ki1 + kj hi) + ?k;] 
= (1—60)*ki. + (1— 00[2k; 1+ k; hi + (ri — 1)kii] 
(1 — 8)?8? [ki + (ri ~ 1) (2k; + ki_ihi) + kia] 
(1 — 0)0° ((r; — 1)k; + 2ki_1 + kj 4 h;] + 05k. 
故 (17) 式 变 为 


Li(z)=(1 — 0)*(fi1 — kii) + (1 — 08A; + (1 — 89?0? B; + (1 — 0)0?C; + 0*5(f, — ki), 


TO — — — — I BR & 学 学 W m 
这 里 
A; = (1 + ri)(fi-i — kii) + hi(di_1 — kj 4), 
B; = ri(fi-i + fi) + hi|di-i — di — (ri — 1)k, 4] — ki — (2ri — 1)ki—1, 
C; = (ri + 1)fi — hi(di + kj 4) — (ri — 1). 
如 果 4; <0, B; «0, C; € Oo, NIE U(t) < 0 对 所 有 ze [rii zi 均 成 立 ， 于 是 该 定理 成 立 。 
定理 4.5 给 出 数据 {(zi, fa di kf kP),4 = 01 onb Afi > kt, ， 则 由 (1) 式 定义 的 
三 次 有 理 插值 样 条 函数 R(z) 在 [zi-1,7i] 上 位 于 二 次 曲线 段 &*(z) 之 上 的 充分 条 件 是 : S 
数 ri(ri > 3) 满足 如 下 不 等 式 组 
(4 ri)(fi-1 — Kia) + (disi — &i*1) 2 0, 
ri(fi-1 十 f) + hi [dii — d; 一 (ri 一 1)k i] 一 k; — (2r; — l)k; , 20, (18) 
(ri -- 1)fi — hildi + k* 4) — (ri — 1)k; 2 0. 
定理 4.6 ”给 出 数据 {(zi fidis es KD RS KT) E01, n] Ek? < f <ko WHR 
定义 的 三 次 有 理 插值 样 条 函数 R(x) 在 [xi_1,zi] 上 位 于 二 次 曲线 段 k(z) 之 下 且 位 于 二 次 曲线 
段 k*(z) 之 上 的 充分 条 件 是 ， 参 数 7; (ri > 3) 满足 不 等 式 组 (16), (18). 
5 ”误差 估计 


B R(z) 是 f(z) 在 [zi-1,24] 上 的 由 式 (1) 所 定义 的 有 理 插值 样 条 函数 ， 当 f(z) < C? a, D] i, 
利用 Peano-Kernel 定理， 可 得 


RIN = f) RE) - [^ FORE- 04er. 
其 中 


1—8)8? [ri (z; —7) - h;]d-0?. i—-T) 
[(z —7)+] | (c-r) CO, meer 
T—T7)4 


一 2 ; nme. n 站 3 ;- 
. (1-8)0 [ri(zi T) hi]+0 Ti D. r«T«Z 


(1-8)2-(1-0)6(r, —1)4-02 


q(T), E<TX< Ti 


H | p(T), Zi-1<T<Z 
考察 作为 变量 7 的 函数 的 核 函 数 Rolle- 7)+] 的 性 质 。 


首先 考察 g(7) 在 [z, zi] 的 性 质 。g(7) 能 被 写成 


02[(zi — x) — (zi — T)((1 — 0)ri+O)) 
(0--G-9vü-SAn-De7 


因为 
(1 — 6)(z; — z)(1— ri) 
(1—89 ro —1) 0? ^ 9 
TN 8?(1 — Os 
q)- —8y1-66(0,-i1):9 ^? 


q(z) = 
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WE g(7) TE [r, 2] 上 的 零点 六 满足 如 下 关系 
(1— 6)h; 
A AEN, 
以 上 情况 说 明 ， 当 r* <r < zi 时 ，gq(7) > 0; mer«r B, 
再 考察 p(7) 在 [zi-tz] 的 性 质 ， 因为 p(7) = (x — 7) + gq(7), 


Zi_1<Z<T* 一 0i 一 < Ti. 


p(z) = q(z) «0, p(zi-1) = (z ~ zi-1) + a(zi1) = Mr 


于 是 ，p(7) 在 (zii, 2) 的 零点 为 
[1 十 6(ri ana 2)]hi 
i 1-60(—-1) 


以 上 情况 说 明 ， 当 zi_1 <T< 办 时 ，p(7) > 0; Ar «r«zH, p(r) «0. XXFE, A 
nun < Wol[ f^ sare f potere f tan [am]. 


- uei [^ omar foie [ aire f^ e] 


* 


hå | eo ris 0) 
21lf ls a 
其 中 
w(r;, 0) = 8?(1 — 6) [£3 (r; — 1)? (ri — 3) — 26?(r? — 3r? — br; +7) 
+ 6(r? — 2r? — 16ri + 15) +r? + 3r; — 6], 
v(ri,0) = [ri - €(r; - 1)] [1 + &(r; - 2] [a — 6? - e(1 — 6)(r; — 1) + €]. 
4 
ax. Ure 9) 
= occi v(ri,0)' 
则 有 如 下 定理 。 


定理 5.1 dE f(x) € C?(a,b, A : a= zo < Ti € c € In = 65 是 区 间 [a,9] 上 的 一 个 分 
划 ，R(z) 是 由 式 (1) 所 定义 的 有 理 插值 样 条 函数 ， 当 z e [enz] A 


| atr < ËI 


6 ”数值 例子 
节 插 值 点 的 导数 值 di 按 单调 数组 和 凸 数组 分 别 采用 误差 阶 均 为 O(h) 的 逼近 方法 ， 
1) 几何 平均 值 通 近 法 
do = AT thha A hh 
di = AP Chia RDUM EUN i—1,2,---,n- 1, 
dc ERA IER. 


"» 00000 I RW X * 学 B — 
2) TORCEXHEOEXT E: 


do = Ai(l + hi/h2) — Ago(hi/ha), 
di 一 (Aihi+ı + hiAigi)/ (hii T hi), i= l, 2, ta,0 — 1, 
dn = Anll 十 hn/hn-1) 5 An,n-2(hs / hi 1), 


其 中 
Anc ELA, Avus fL hoch 


Qs 


X2 — To Tn 一 Zn—2 
例 1 f(z)-l,ze[-2,-02], HEARA z = [一 2, 一 1.5, —1.0, —0.8, —0.6, —0.4, —0.2]; 
导数 值 按 几 何平 均值 逼近 法 所 得 图 像 为 图 1， 其 中 ， 左 图 取 w = 10, HRR qi = 24。 从 图 1 中 
可 以 看 出 ， 只 要 满足 (10) 式 ， 就 可 得 到 保单 调 插值 。 


图 1: 保单 调 有 理 样 条 插值 


例 2 数据 点 
z = [一 2, —1.5, —1, —0.5,0, 0.5, 1, 1.5, 2, 2-5], 
y = |-8, —3.375, —1, —0.125, 0, 0.125, 1, 3.375, 8, 15.625]. 


图 2 中 ， 左 图 取 gq; = 2， 右 图 取 qi = 15， 是 采用 导数 的 算术 平均 值 逼近 法 ， 给 出 的 有 理 样 条 插 
值 图 形 ， 调 整 g; 满足 (11) 式 即 可 得 到 保 凸 插值 。 

例 3 设 f(z) = cos2( 舍 )。 取 插值 结 点 0.0, 0.75, 1.5,2.25,3.0, 3.75,4.5,5.25,6.0， 即 等 距 划 
4j, h = 0.75， 其 在 [0,6] 上 R(z) 为 由 (1) 所 定义 的 有 理 样 条 曲线 。 取 上 约束 曲线 g(z) 和 下 约 
东 曲 线 g*(z) 分 别 为 


11-2z, 0<z<1.5, 09—2z, 0<z<1.5, 

3z—0.9, 15<z<3, ， 3r—11, 15<7<3, 
g(z) = g'(z) = 

3.1 — 2z, 3£€z«4.5, 2.9 一 27, 3<Z<4.5， 


#7—2.9, 45<z<6 2r—3., 45 € r X6, 
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则 当 在 第 一 个 子 区 间 [0,0.75] ERr = 4， 在 其 它 子 区 间 取 7; = 6 时 ， 插 值 曲 线 R(x) 被 约束 于 
折线 g(z) 与 9*(z) 之 间 。 图 3 表示 了 插值 曲线 Ra) 完全 被 约束 于 折线 g(z) 与 9*(z) 之 间 。 


图 3: £55 R(z) 于 两 折线 之 间 
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The Approximation Properties of a Cubic Rational 
Spline Interpolation 


LIU Lin, TANG Yue-hong 


(Department of Mathematics, Nanjing University of Aeronautics & Astronautics, Nanjing 210016) 


Abstract: In the geometric shape design, the shape preserving interpolation of a curve is an important 
topic. A new method for rational spline interpolation of order (3, 2)* (k = 1,2) with control parameter 
and prescribed poles is constructed. For given data satisfying monotonicity or convexity, a shape 
preserving interpolation function can be generated by choosing suitable parameters. The issues on 
shape control are discussed, and the sufficient conditions that constrain the interpolating curves to 
be above, below or between the given broken lines or piecewise quadric curves are derived. With the 
Peano-Kernel theory, the approximation properties of this rational spline are studied. The numerical 
examples show the effect of these methods. 

Keywords: rational interpolation; spline; constraint interpolation; approximation 
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